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Abstract. We consider the Gross-Pctacvskii equation in 1 space dimension 
with a n-well trapping potential. We prove, in the semiclassical limit, that 
the finite dimensional eigenspace associated to the lowest n eigenvalues of 
the linear operator is slightly deformed by the nonlinear term into an almost 
invariant manifold Ai. Precisely, one has that solutions starting on M, or 
close to it, will remain close to M. for times exponentially long with the inverse 
of the size of the nonlinearity. As heuristically expected the effective equation 
on M is a perturbation of a discrete nonlinear Schrodingcr equation. We 
deduce that when the size of the nonlinearity is large enough then tunneling 
among the wells essentially disappears: that is for almost all solutions starting 
close to M their restriction to each of the wells has norm approximatively 
constant over the considered time scale. In the particular case of a double well 
potential we give a more precise result showing persistence or destruction of the 
beating motions over exponentially long times. The proof is based on canonical 
perturbation theory; surprisingly enough, due to the Gauge invariance of the 
system, no non-resonance condition is required. 



1. Introduction 

In this paper we study the dynamics of low energy states of the one-dimensional 
Gross-Petaevskii equation (hereafter also called nonlinear Schrodinger equation, 
NLS) 

1 ^(i)| t=0 = ^(i)eL 2 (i), U°\\^ = i, 1 ] 

where a is a positive integer number and 

d 2 

H = -h 2 -i^ + V, x e R, (2) 

is the linear Hamiltonian operator and V(x) a n-well potential. By this we mean 
that V has n nondegenerate distinct minima x\,...,x n where the potential has 
essentially the same behavior (e.g. one can assume that its first r derivatives are 
equal at all the minima, for some positive integer r > 4). We also assume that 
the potential is trapping, i.e. V tends to infinity as \x\ — > oo. The situation 
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Figure 1. Plot of a trapping potential with n wells. 

we have in mind is that of a Bose Einstein condensate trapped by an unbounded 
potential and also subjected to a periodic- like force field (see Figure^; in such a 
case the parameter e can be thought as a measure of the number of particles in the 
condensate (see, e.g. USD- 
Consider first the linearized problem. The fundamental state of Hq is approxi- 
matively degenerate in the sense that, denoting by 

Ai < A 2 < ... < A„ < ... 

the eigenvalues of Hq, one has 

A n — Ai <c A n _|_i — A„ 

in the semiclassical limit, i.e. fi«l. Then the most interesting situation occurs 
when the normalized cigcnfunctions ip±, ...,ip n , corresponding to Ai, A n , are delo- 
calized among the wells. Indeed, in such Ob CclSG ci solution in <E>o :=span(y>i, tp n ) 
performs a quasiperiodic motion and the probability of finding the particle in any 
fixed well undergoes great changes over a time scale of the order of T = ttH/lu 
with u> = (A„ — Ai)/2. In the case of double well potential such a phenomenon is 
usually known as beating motion and the beating period is given by T := Tth/u), 
with lu = {X-2 — Ai)/2. 

The main question is the behavior of the system when the nonlinearity is restored. 
In the case of a double well potential the problem was tackled in a series of papers 
dHZHiniiniiniEniEIlElESlEZl; in particular, it was shown that, up to times of 
order T, the dynamics is well described by an Hamiltonian integrable system with 
two degrees of freedom obtained by restricting the Hamiltonian, i.e. the energy of 
the system, to $o- In particular, this result has been used in order to show that the 
beating motion is generic for values of the nonlinearity strength e below a certain 
threshold value, while new localized states appear for larger nonlinearity strength 
(i.e. as the number of particles of a Bose-Einstein condensate increase) and for even 
larger values of the nonlinearity strength the beating motion disappears. 
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In the case of a multiple well potential the situation was studied e.g. in the 
paper J2SJ where the authors deduced (formally) the discrete nonlinear Schrodinger 
as an effective equation for the dynamics in <&o and used it in order to study some 
of the features of the model. We are not aware of rigorous results in the case of a 
multiple well potential. 

In the present paper we study the nonlinear dynamics using the methods of 
Hamiltonian perturbation theory for PDEs B] HI G31 and we prove that the 
manifold <E>o, which is invariant for the linear dynamics, is only slightly deformed 
by the nonlinearity into a new manifold M which is approximatively invariant for 
the complete dynamics. By this statement we mean that solutions starting on M, 
or close to it, will remain close to M for times which are exponentially long with 
T; this result improves the one given by [3121312] for double well potentials where 
the control of the approximate solution was given for times of order T. Moreover 
we show (see LemmaEl that the dynamics on M. is described at first order by a 
discrete nonlinear Schrodinger equation 

#j = + + A j-i^j-i + #j l^jl" , 3 = l,-,n, (3) 

where ipo — tpn+i '■= 0, Aj, Sj,rj are suitable constants and where, having in mind 
the case of Bose-Einstein condensates, | 2 represents the fraction of particles in the 
j-th well for j = 1, n (see Sect. I2.4l for their precise definition). In particular, it is 
quite easy to study the discrete NLS © from the anticontinuum limit ^1] obtaining 
that when r\ is large enough then solutions corresponding to almost all initial data 
have the property that \ipj\ 2 is essentially a constant of motion. In the case of 
double well potential, we get that the dynamics on M is (up to an exponentially 
small error) that of an integrable Hamiltonian system with two degrees of freedom. 
This allows us to control also the trajectories of the solutions on M. showing that 
for small r\ the beating phenomenon persists for exponentially long times, while for 
higher values of r\ only motions which are essentially localized in one of the two 
wells exist, at least for exponentially long times. 

We emphasize that from the technical point of view the main result (Theorem|5J) 
is quite surprising since the application of canonical perturbation theory is typically 
possible only when some non-resonance conditions are satisfied. On the contrary, 
here the result is valid for any multiple well potential, whose eigenvalues might 
fulfill arbitrary resonance conditions (i.e. the eigenvalues can be linearly dependent 
over the relative integers). This is possible since NLS is an infinite dimensional 
Gauge invariant Hamiltonian system. To explain how this property is exploited we 
recall that canonical perturbation theory allows us to remove from the Hamiltonian 
all non-resonant monomials. In particular, given an arbitrary monomial it can be 
eliminated if it is non-resonant. However, in NLS only Gauge invariant monomials 
appear, and we will show that, for a Gauge invariant monomial, the non-resonance 
condition is (almost) trivially fulfilled. Actually, in order to avoid any restriction 
on the potential we have to use a resonant construction which is not self evident a 
priori (see Subsection 13. 2(1 . We think that these ideas could be useful in the study 
of Hamiltonian systems with symmetry, and possibly also for the investigation 
of further dynamical properties of NLS. A further technical ingredient which is 
fundamental for the proof is the use of Sobolev like spaces constructed as the 
domains of the powers of Hq . To use such spaces one has to show that they form 
Banach algebras under the pointwise multiplication. Here we give a detailed proof 
of such a property that we think could be useful in further investigations of NLS. 
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The paper is organized as follows. In Section [3 we state our main results (The- 
orem Theorem |3 and their corollaries). Section [3 is divided into 5 subsections: 
in Subsection 12 . II we review some known facts about the structure of the low lying 
eigenvalues of the linear Schrodinger equation with an n well potential; in Subsec- 
tion l2.2l we introduce the Sobolev like spaces in which NLS equation will be studied 
and give their main properties; in Subsection 12 . 31 we will state the result on the ap- 
proximate invariant manifold; in Subsection 12.41 we will give the effective equation 
on the approximatively invariant manifold and deduce the localization properties of 
the solutions; finally in Subsection 12 . 51 we will study the particular case of a double 
well potential. In Section [3] we prove our main results. This section is also divided 
into 6 subsections that correspond to the different parts of the proof. The proof of 
Theorem 1 (algebra property of Sobolev like spaces) and of some technical Lemmas 
are left to appendixes A and B respectively. 

Acknowledgments. DB would like to thank Panos Kevrekidis for pointing to his 
attention the paper [231 an d the connection between the Gross-Pitaevskii and the 
discrete NLS equations. 

2. Main results 

2.1. Linear theory. 

Hypothesis 1. The potential V(x) € C°°(R) is a real valued function such that: 

i. V(x) admits n minima at x\ < X2 < ■■■ < x n such that 

V{x) > Vmin = V(xj) = 1, V.x el, i/ Xj, j = 1, ...,n; (4) 

ii. There exists a constant C > 1 such that 



C-\x) 2 < V{x), (x) = y/l + x 2 ; 
iii. There exists a positive m > 2 such that for any k £ 
d k V{x) 



dx k 



< c k ( x y 



for some positive constant Ck; 

iv. The minima are nondegenerate and 

^ = C>0 (5) 

with C independent of j; 

v. The shape of the potential at the bottom of the minimum Xj is approxima- 
tively independent of j; precisely: there exists r > 4 such that 

d k V d k V 

Hereafter, C will always denote a positive, typically large, constant whose value 
changes from line to line, and which is independent of h, e and t. 

The operator Hq formally defined by @ admits a self-adjoint realization (still 
denoted by Hq) on L 2 (R) (Theorem III. 1.1 in |S]) with purely discrete spectrum. 
Let Afe, k £ N, be the non degenerate eigenvalues of Hq 

Ai < A 2 < A 3 < A 4 < . . . < A fe < . . . , X k oo, 

with associated normalized (in L 2 ) eigenvectors ipk(x); the set {(fikix)}^^ is an 
orthonormal base of L 2 . 
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The lowest part of the spectrum can be studied in the semiclassical limit using 
the construction of ^2], that we shortly recall. 

Having fixed a positive constant a > 1 = V m i n we consider the set V~ 1 ((— oo, a)) 
and we assume that a is such that this set is the union of n disjoint open sets Uj 
with Xj G Uj . Having fixed j we consider the operator formally defined on L 2 as 

,o d 2 



H, 



dx 2 



where Vj is a modified potential defined as 



V, 



for x S 



V{x) 

max[a, U(ir)] for x Uj 



(7) 



Let Xj be the lowest eigenvalue of Hj with associated normalized eigenvector ipj. 
Using the semiclassical construction of the eigenvalues close to the bottom of a well 
(see e.g. 01221) one has that assumption Hyp. [T\v., implies that 



and 



for any i,j = 1,2, 
Let 



Xj =1 + 0(H), \Xj-Xi\ < Ch r/2 



\\<fij{x) - (pi(x + x t - Xj)\\ L 2 < Ch r/2 , 



x j + 1 



y/V(x) - ldx, 3 = 1,2,. 



1. 



be the Agmon distance among the two minima Xj and ij+i. Let T any fixed 
positive real number such that T < mhij Tj . Then, from the theory of |12| there 
exist some constants cy such that for any i and j 



22,C k j<pk 
k=l 



< Ce- T ' h , 



\\vm\\ L ~<Ce- T ' h , i + 3 , 

\Kj — A jl < Ce~ v l n for some ij G {1, ...,n} . 



(8) 

(9) 
(10) 



Moreover, the projector n on $0 :=span((^j) is a bijection among $0 :=span(iy9j) 
and $0 itself. Remark also that one can choose the functions (pj to be real valued. 
Hence, the lowest n eigenvalues of Hq fulfill 



C^H < Xj : - 1< Ch , 3 = 1, 



for some C > 1 and 



= \{X n - Ai) < CK-' 2 



Furthermore, making use of the same arguments, it follows also that 



inf 



inf 



j=l,...,n Ae<r(J? )-{Ai,...,A n } 

We also use the following projectors 



[A - Xj] > c- x h 



(11) 
(12) 



n 



(ipj,-)ipj and II C = I — n and $ = ILL 2 ; 
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Remark 1. As already emphasized in the introduction the most interesting situ- 
ation occurs when the true eigen] 'unctions (fj are delocalized between the wells, a 
property that occurs if the potential is exactly periodic in some region or more gen- 
erally if the order of magnitude of some the off diagonal elements of the matrix 
formed by the constants Cij is of the same order of magnitude as the difference 
between the approximate eigenvalues Xj , a property that is guite difficult to ensure 
in a general situation. For this reason, essentially in order to fix ideas, we decided 
to state our assumptions in the form (i-v) above. 

2.2. The nonlinear system: Analytic framework and well posedness. 

Definition 1. For any integer s > define the Hilbert space X s := D(H S J 2 ) 
endowed by the graph norm; more precisely in X s we will use the following norm 
equivalent to the graph norm 

H\ 

The main step for the proof that the spaces X s form a Banach algebra under 
the pointwise multiplication is the following Theorem. 

Theorem 1. Let s be any positive integer number. For small enough h the two 
norms 

Ml and \\(-h"Ay/^\\l 2 + \\V^m 2 (14) 

are equivalent with an h independent constant. 

The proof, which is a semiclassical variant of the proof of Lemma 7.2 of is 
deferred to Appendix lAl 

Remark 2. In particular one has that a function <f> is in X s if and only if it belongs 
to the Sobolev space H s and it decays at infinity so fast that \<fi\ 2 V s is integrable. 

In the spaces X s with s > 1 the system Q is semilinear, since, using (|14|) and 
Gagliardo-Niremberg inequality one has 

Corollary 1. For any integer s > 1 there exists a positive constant C s independent 
of h, such that 

HifaWs <C.*~ 1/2 II&IUI&L> V^^er,k<i. (15) 

Moreover, the map 

X s x X s 3 (4>,4>) i-> |0| 2<T </> S X s 
is entire analytic map for any h > and fulfills 

||l0l 2 ^|| s <a^il0llf +1 - (is) 

Remark 3. In dimension d > 1 this result remains valid provided s > d/2. 

Then by standard Segal theory (see e.g. ^B]) the system JQ) is locally well posed 
in all the spaces X s with s > 1. Actually s > d/2 is enough, but our proof only 
applies to integer values of s; in jJBJ a Strichartz inequality argument was used to 
show that it is also locally well posed (LWP) in X s with some s smaller than d/2. 

From now on we assume that the index s of the space is a fixed positive integer 
number and fulfills the condition s > 1. In the following, in order to fix ideas, one 
can just think of the case s = 1. We will denote by d{.\ .) the distance in the norm 
of X s . 



H, 



L 2 



' L- 



G X s . 



(13) 
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2.3. The nonlinear system: Approximatively invariant manifold. In order 
to state our main result we assume that the size of the nonlinearity is small enough, 
i.e. |e| <C h? , and we introduce the small parameter. 



(17) 



where to was defined in 

Theorem 2. Consider the system Q and fix a positive s > 1. There exists a 
positive /i* such that, if [i < fj,*H 3 / 2 , then there exists a manifold A4 (dependent on 
all the parameters of the system) with the following properties: 

i. M. is close to <&0; *-e. 



d($ ,M) < C 



(18) 



where 



n 



d($ , M) = sup inf | 

and where \\(f\\ s is the norm (|13H . 
ii. Let 



d = d(i>°,M) = inf H°-<p\ 

<p£M 



be the initial distance from M. and let 



5 = max ^ do , exp 
Then for all times t fulfilling 



[i 



,fi 3 / 2 



2// 



t < 



ChuS 



one has 



d(ip\M) < C5 



< c 



jf_ 

ft 3 / 2 



(19) 

(20) 

(21) 
(22) 



Such a manifold M is called an approximatively invariant manifold. 
Remark 4. The most interesting cases are when 



5 = exp 



2/x 



S = C 



h 3/2- 



Indeed in the first case the time of validity of all the estimates is exponentially long, 
while in the second case it is easy to obtain the following Corollary. 

Corollary 2. Assume that ||n c ^°|| s < C/ifi 3 / 2 then, up to the times 

h 3/2 



t < 



(23) 



the estimate 1221) holds. 



s 
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Remark 5. This Corollary is a direct extension of the results of [111 1201 121) in 

which the estimate ()22ll has been proved (for the double well potential) for a time 
scale of order 

T = — < C— -C ^ 2 h 3 ' 2 . 

LU [1 

The improvment is due to the fact that our construction implies that M. is linearly 
stable up to an exponentially small error. 

2.4. The nonlinear system: discrete NLS and suppression of tunneling. 

To start with we remark that NLS is a Hamiltonian system (see Subsection 13 . II for 
a precise description) with Hamiltonian function given by 

£(V,$:=£o(^)+£Po(^$ (24) 

where 

£b(t/>,$) := / tp(x)H tp(x)dx 
Jr 

and 

:= -^-r\\r +1 \\h ■= — I r +1 r +1 dx 

a + 1 - ff + lJi 

The general idea is that for initial data close to $o the system should be well 
described by the Hamiltonian system obtained by restricting £ to $o which is close 
to $o- Denote 

n 

V> = E^ e $0 
3=1 

then we have the following 

Lemma 1. T7ie restriction of (I24|l to $0 tofces i/ie /orm 

n n n 

% = E( n + ^)i^i a + ec E i^-i 2ff+2 + E ^ fe^-i + ^^--i] + 

3=1 i=i 3=2 

+ 0(e- 2r / ri ) + eO( e - r / ri ) 

where T was introduced before equation (JSJ and 

1 " /■ 
Q = - V X 3 , vj := / <pjH <Pjdx - = 0(?T /2 ) 

71 3 =1 •'K 



3=1 



Cj := / y>j 



jHo^dx = 0{e- T ' h ) 
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Proof. Indeed, we have that 



^ l*o = f ol|, +e^o| 4o - 
The first term takes the form 

n 

£ °ko = E \^i\ 2d n + E + E ^i d v 

3=1 |i-j|=l I<-J'|>1 

with 

dij = / (fitHoipjdx 
Jr 

and where (see, e.g., Q2HHI) = 0(e _r/fi ) when \i-j\ = 1 and d hj = 0(e- 2T l h ) 
when \i — j\ > 1. For what concerns the second term, following and making 
use of ©, we have that 

n 

° 3 = 1 

where = O (fi.~ CT / 2 ) is approximatively independent of j. Remark also 

that, since the functions ipj are real valued, the quantities Cj turn out to be real. □ 

Hence, up to higher order terms, to the Gauge transformation ip — ► e int l h ij) and 
to a rescaling of time t — > ojt/H the restriction of l|24|) to <&o is given by 

n n n+1 

/Co :=X>K ; j! 2 + l^| 2CT+2 + E A i(^j-i + ^i-i). (25) 

3=1 3=1 3 = 1 

where 

?7 = — = f-feL) , *7 := — = 0(1) , A,- := ^ = 

uj \ujh' J i z J uj uj \ uj 

and the equation of motion of (|25J) are given by 

w, = — i — =-, with 7 = 1, n. 

V3 diPj' J 

The system (|25|) has an integral of motion (the restriction of the square of the L 2 
norm to <E>o) given by 

n 

I:=El^f (26) 

3=1 

We analyze now the consequences of our Theorem for the dynamics. From the 
proof of Theorem[21we will be able to exactly describe the restriction of the system 
to M. . Actually, we can state the following Theorem^ which is a result of Theorem 
Instated below. 

Theorem 3. Under the same assumptions of Theorem^ there exists an analytic 
canonical transformation T :IA — ► X s , with hi an open neighborhood o/$o with size 
independent of fi and h, such that 

£oT = m + ojJC+n 

where 
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i. JC o II = JC, i.e. IC depends only on the variables tpi, ip n , 

ii. {T, IC} e0 i.e. I is an integral of motion for the system with Hamiltonian 
JC 

iii. JC = JC Q + 0{p/h z ' 2 ) 



1Z= 



exp — (y^j^ — ^ + ( || n c ^ || ^ and similar estimates hold for its 
vector field. 

v. T = I + 0(/x/fi 3 / 2 ) i.e. the transformation is close to identity. 



Remark 6. In this framework the manifold M. turns out to simply be A4 = T($o)- 

In particular it follows that the dynamics on M. is, up to a small error, the same 
of a Hamiltonian system with Hamiltonian function close to JCq with an integral 
of motion given by I. Thus, it is possible to deduce the following Corollary of 
Theorem below, which is particularly relevant in the double well case. 

Corollary 3. Under the same assumptions of Theorem\^one also has 



\l(t)-l(0)\<C 
up to the times (|2L)[1 . 



£3/2 



\K (t)-K (0)\<C 



h 3/2 



(27) 



We focus now on the exponentially long time scale, thus we assume that the 
quantity 5 of Theorem [5] is given by 



5 = exp 



2^ 



and from (|20[1 it follows that up to times of order e 
fulfills the a priori estimate 



the vector field of IZ 



1**1 



< C exp 



/'■* 



Thus, up to an exponentially small drift, a Gauge transformation and a rescaling 
of time the dynamics on M = T(<i>o) is that of the n-dimensional Hamiltonian 
system JC, which is a perturbation of JCq, i.e. of the discrete NLS l|25l) . 

It is worth mentioning that when ?; = the dynamics of JCq is that of n decoupled 
harmonic oscillators corresponding to the normal modes of the linearized system. 
If Aj <§C Si, for any i and j, then the normal modes are localized, i.e. each normal 
mode essentially involves only one of the ipj; on the contrary, in the much more 
interesting case where the Sj are of the same order of magnitude of the A^, typically 
the normal modes are collective motions of the system and, correspondingly, in 
typical solutions the term |^j(f)| 2 undergoes great changes for each j. 

In the opposite limit rj — > oo (anticontinuum limit, see |14j1. JCq becomes a system 
of decoupled anharmonic oscillators. Correspondingly | V^j? 1 2 is a constant of 
motion. One can use KAM or Nekhoroshev theory in order to study the dynamics 
of JC when 1 rj < oo and to deduce results on the dynamics of the complete NLS 
equation. 

Here we will state a result that can be obtained in this way. To this end, for 
any p > fixed, we will denote 



s!: 



3=1 



= 1 and |^| > p , Vj 
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and we will denote by <S"| its Lebesgue measure. 



Theorem 4. (KAM theorem IC) Consider the Hamiltonian system IC, then there 
exists a constant rj^, such that, for any \i]\ > there exists a set S v C 5™ with 
Lebesgue measure estimated by 

\S V \ > \S n p \-C V - 1/2 (28) 
with the property that, if the initial datum is in S-q then the solution of the Hamil- 
tonian system with Hamiltonian function IC is quasiperiodic and fulfills 



< Cn 



-1/2 



(29) 



for any t. 

To state a corresponding result for the NLS equation, consider the set C of the 
ip G X 3 having L 2 norm equal to 1 and fulfilling 



d(ip,M) < Cexp 



2[i 



For ip G C, denote Y^j=i ^iVj — ^P an d 

|^-| > p , Vj} . (30) 
Corollary 4. Consider the NLS equation 1(1)1 . under the same assumptions of 



s n p ~ e c 



theorem \^ assume also r\ large enough; then there exists a set S n C <S" whose 



"measure" is estimated by 



> 



us; 1 



-1/2 



such that, if ip £ S v then along the corresponding solution one has 



^■(*) - Mo) 



-1/2 



(31) 



(32) 



for all the times t fulfilling l|2U|l . 

Remark 7. All the constants in the above Theorem^ depend 
can find some n independent statements. 



in |21 El one 



2.5. The double well potential. In the particular case of a double well potential 
one can get much more precise results, both for the linear and for the nonlinear 
system. 

By a double well potential we will mean here a potential V(x) E C°°(R) fulfilling 
assumptions (i-v) of Section ^. II (with n = 2) and which moreover is symmetric with 
respect to spacial reflection V{— x) = V(x). It is well known that the splitting u> 
between the two lowest eigenvalues fulfills the asymptotic estimate 

u= l -{\2-\i)<Ce- T l\ (33) 

for any T < Ti where Ti is the Agmon distance between the two wells and C is a 
positive constant (depending on T), and moreover the normalized eigenvectors <$\,i 
associated to A^a can be chosen to be real-valued functions such that ip\ and if2 
are respectively of even and odd-parity, thus, defining the single well states 



PR 



<p 2 ] and 



<PL = ^ to 



V2] 
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they essentially coincide with the functions <pj used in the multiple well case. In 
particular they fulfill 

HVBVilU- <Ce- T ' h , T>0. (34) 
Thus, for ip S $o, in this section we will write 

ip = iPiVl + 1p2fR (35) 
Here the tunneling gives rise to the so called phenomenon of beating: for almost 
any initial datum ip € <E>o the expectation value of the position 

{x}* = [ x\^(x)\ 2 dx (36) 

periodically oscillates between positive and negative values with a period given by 
T := ttH/lu. 

In this case the restricted approximate Hamiltonian /Co takes the form |21| 

/Co := faih + V^i + v (l^i | 2 " +2 + \H 2a+2 ) (37) 
with the same definition of T] as in the previous subsections. Then Theorems [21 
and hold together with their corollaries. The main improvement that one can 
obtain in the double well case are due to the fact that since the system /C is now a 
system with two degrees of freedom with an integral of motion independent of the 
Hamiltonian (namely I), then it is integrable. This allows us to describe in a very 
precise way the trajectories of the system K which are just the intersection of the 
level surfaces of the functions /C and I. This is possible since JC is close to /Co- 

To be definite, from now on, we will restrict to the case a — 2. 

The system /Co, cf. ((23, has been already studied in [TT] (see also [T71 I2H 
obtaining that, for \rj\ < 2 almost all solutions perform beating motions, while 
at t] — ±2 a bifurcation occurs and new equilibria, localized close to the minima 
of the Hamiltonian function, appear. As r\ increase the domain of stability of 
such solutions increase its size, so that, for r\ large enough essentially only localize 
motions exist. Concerning the complete system we can state that if r\ is not at a 
bifurcation point, then non-homoclinic trajectories associated to the Hamiltonian 
/Co approximate the solution ip* for times of the order (|20l) . 

Corollary 5. Under the same assumptions as in Theorem\^ assume also i] ^ ±2; 
consider an initial datum such that ipi(0),tp2(0) ^ 0, and 

MO) + \ - 1 if v>o 
Mo) i - I if v<o 

and also such that S < C/ih 3 ^ 2 ; then, there exists a positive constant depending 
only on how much the above quantities differ from the considered values, such that, 
provided \x < [i^h?!" 1 , there exists a solution of the Hamiltonian system (|25|l with 
trajectory 7 such that 

for the times (|23Jl . 

Remark 8. Homoclinic trajectories are absent when \rj\ < 2; initial conditions such 
that /Co(0) = § — 1. for r\ > 2, and /C(0) = 1 ~ §' f or V < ~^ corresponds to an 
homoclinic trajectory. 
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Remark 9. The topology of the trajectory 7 is determined by the condition 

d{^^)<C^ (39) 

in the sense that all curves fulfilling this condition have the same topology if the 
assumptions of the Corollary are fulfilled. 

Remark 10. Thus one has that also for the true system beats are present for 
\rj\ < 2 while their importance decreases as \n\ increase above 2. In particular for 
large values of rj only motions localized close to one well are present |25| at least 
for the time scales controlled by our theorems. 



3. Proof of the main results 

In order to simplify the notations all the proofs will be carried out in the case of 
a potential with only 2 wells. 

3.1. Hamiltonian Formalism. First consider the real Hilbert space 

x r : = -D((#o,r) s ) 

where i?o,K is the operator Hq restricted to real valued functions. 

We make X^ © X^ a symplectic space by introducing the semiclassical symplectic 
form 

a((p,q)\(p',q')) ■= H [p'{x)q{x) - p[x)q'{x)] dx. 



Given a smooth real valued function TL(p, q), then we define its Hamiltonian vector 
field Xu € X£ © X£ by the property 

a(X H ,h)=dTLh (40) 

for any h = (h p , h g ) € X£ © X£, where dTt denotes the differential of TL. It is well 
known that Xu is in general defined only on a subset of X£ © X£ . Define also the 
L 2 gradient W P TL of TL with respect to p by 

(V p H,h p ) L2 := f \7 p H{x)h p (x)dx = d p Hh p , Wh p € (41) 

JR 

and similarly we introduce the quantity V q TL. Then 

x n = n- 1 (-v q n,v p n), 

and thus the Hamilton equations of TL are given by 

j f (p, q) = X n (p, q) (p= -\^ q H , q = ]^ p h) . (42) 

The Poisson brackets between two functions TL and K, is defined as 

{TL,IC} :=a(X n ,X,c) = -\ [ [V P TLV q IC - V q TLV P IC} dx (43) 

" JR 

which in general is only defined on a subdomain of X£ © X£. 

We shall use complex coordinates in XjZ © X£ identifying this space with X s , 
through 

(p,q) ^>ip= -j=(q + ip). 
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Therefore, we set 

V^ = -J|(V g -iV p ) and V v - = -L(V 9 + zV p ) (44) 

so that, if "H = T~C(tp, tp) is a smooth real valued function, we have the identification 

Xnty,® = ~VjH(i/>,<$) . (45) 
and in complex coordinates the Poisson brackets are computed by 

{H, K.} := \ [ [V^V^/C - V^V^/C] dx (46) 

With such a notation then the NLS QJ can be written in the form of a Hamil- 
tonian system, the corresponding Hamiltonian function being the energy eq. (|24|) . 

Remark 11. Such a Hamiltonian is invariant under the action of the Gauge group 

4>{x) -> ^(^)e~ J/5 , (47) 

for any (3 independent of x. The corresponding conserved quantity is the L 2 norm 

AT(tl>,$) = [ \^P(x)\ 2 dx. (48) 



Equivalently one has 

{S,Af} = 0. (49) 

Let {\k\kL\ and {<Pk( x )}kLi be the eigenvalues and the normalized eigenvectors 
of H , let 

oo 

*p(x)=J2&Mx) , (50) 
fe=i 

and define the Hilbert spaces t\ of the complex sequences such that 

IIC||*:=$>£|&| 2 <00 ( 51 ) 

k>l 

In such a way we have defined the correspondence 

1>€X'~Uty) :=c=(Ci,C2...0,-)e^ 

which is a unitary isomorphism. In particular, if £ is the Hamiltonian 1)24)). then 
£ oW _1 (still denoted by £) is the Hamiltonian of the same system written in terms 
of the new variables £• In terms of these variables the quadratic part £q of the 
Hamiltonian is given by 

£> = X>lO| a ( 52 ) 
and the Poisson brackets can be written as 



fc=i 



90s 9Cfc aCfc 5Cfc 
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Remark 12. From now on we will work in the space £ 2 and moreover, in order 
to simplify, we rescale time by the transformation 1 1— > t/h, thus the Poisson 
brackets take the form 



" r dH dIC dU dIC 



k=l 



(53) 



3.2. Non coupling monomial. It is useful to introduce also a different notation 
for the first two variables (here we recall that we are working in the double well 
case) and for the remaining ones, thus let us denote 

til := Ci . U2 ■= C2 , Zj := , j > 3. (54) 
Consider now a monomial of the form 

( K ( L = u k u l z m z n (55) 

where we used the notations 
K = (fc,rn), L = (l,n) 

k = (ki,k 2 ), m = (m 3 ,TO4,TO 5 , ....), I = {h,h), n = (n 3 , n±, n 5 , ....) (56) 



,k — „,fcl„,k' 



i„, K 2 ~m 



Remark 13. A monomial of the form (15511 is Gauge invariant, i.e. invariant under 
the transformation 

ip h-y e^tp that is Cfc ^ e i/3 (jt 
if and only if \K\ — \L\, where 

\K\=Y,\KA- 

3 

In fact (|55|l is Gauge invariant if, and only if, 

= {M(C, 0, C K C L } = * (y^iKj - Lj)j ( K ( L (57) 

where j\f ((, , C) = Si>i | Cj" 1 2 * s ^ e L 2 norm; that is \K\ — \L\ again. 

Due to our assumption Hyp.l on the potential one has 

Lemma 2. Let o~(Hq) be the spectrum of the linear operator Hq and let h be small 
enough. There exists a sequence of (not necessarily continuous with respect to h) 
functions {£' 7 (?i)} 7e N, < Eq < 1 and there exists a positive constant C > 1, 
independent of H and 7, such that: 

i. [E 7 - C~ x h, E 7 + C^h] n o-(Hq) = 0; 

ii. 1 < E 7 — E-y—i < 3 for all h and all 7 > 1; 

iii. For any h fixed, we consider the sets of indexes 

J 1 (K) = J 7 := {j e N : E y -i < Xj < E 7 } , 
then the cardinality of these sets is estimated as 
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Proof. The proof is an immediate consequence of the following result (see Theorem 
(V-ll) in [THj, see also Theorem (XIII-81) in {TQ): 

%« = ^[|{(a;,p)eR 2 : a<p 2 + V(x)</3}\+0(h)] 

where iVr ai #i is the number of eigenvalues of Hq contained in the interval [a, [3] 
and where |.| denotes the Lebesgue measure of a set. Indeed, let us consider the 
intervals [27,27 + 1], 7 £ N, then the number of eigenvalues of Hq belonging to 
these intervals is given by 

N{i) := iV [27 , 27+1] = 

= i^ (l ' p)eR2 : 2 7<P 2 + ^(^)<27+l}|+0^)] <j 

since V(x) > C(x} 2 , for some positive constant C independent of h and 7. From 
this estimate it follows that there exists at least one value E 1 E (27, 27 + 1) and 
C > 1 satisfying L. Furthermore, conditions ii. and iii. immediately follow since 

1 = (2 7 ) - [2( 7 - 1) + 1] < E 1 - £ 7 _i < (27 + 1) - 2( 7 - 1) = 3 

and 

T 

□ 



# J 7 < N{i) + JV( 7 - 1) < — . 



We fix a sequence with such properties. 

Hereafter, all the perturbative construction will involve only Gauge invariant 
monomial (hence \K\ = \L\). Keeping this in mind we give the following 

Definition 2. A monomial of the form (I55|l will be called coupling if the indexes 
(K, L) fulfill the following conditions 

i. \K\ = \L\ 

ii. |m| + \n\ = 1,2 

iii. if mi — 1 and nj — 1 then i £ J 7 and j £ Jy with 7 ^7'. 

A monomial which is not coupling will be called non coupling. A polynomial 
containing only coupling (reps, non coupling) monomial will be called coupling 
(reps, non coupling). 

Remark 14. For Gauge invariant monomials the condition i. is always fulfilled. 
Furthermore, in terms of the indexes k^l^m^n condition i. reads 

\k\ - \l\ = \n\ - \m\ (58) 

Remark 15. Any Gauge invariant analytic function can be uniquely decomposed 
into the sum of a coupling and a noncoupling part. We recall that, as in finite 
dimensional spaces, an analytic function is a function whose Taylor series is con- 
vergent (see e.g. |15| ). 

We also define a new {h dependent) norm in the space X s as follows 

Definition 3. Denote 

M E {z,z):=Y,E s 1 £ M 2 (59) 

7 jeJ T ,j>3 
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then the quantity 

\\(u,z)\\ 2 E = \\Cf E := K| 2 + \u 2 \ 2 +M E (z,z) (60) 

will be called E-norm of Q. Denote £ E the space of sequences z = {zj}j>3 equipped 
with the norm Af E (z,z). By abuse of notation sometime we will also write 

\\ Z \\ E -=Ne{z,z). 

Remark 16. The E-norm is equivalent to the standard norm of X s with an h 
independent constant. The proof is a trivial computation and is left to the reader. 

Lemma 3. Let ( K ( L be a non coupling Gauge invariant monomial of degree at 
most two in z,z (i.e. \n\ + |m| < 2), then one has 

{M E ,Q K Q L }=Q and {\u l \ 2 + \u 2 \ 2 X K C L } -0 (61) 

Proof. First remark that if u k u l z n z m is noncoupling of degree at most two in z, z 
then one has \k\ = \l\ and thus 

{uiui + u 2 u 2 , u k u l z n z m ) = i{\k\ - \l\)u k u l z n z m = 0. (62) 

Furthermore, either it is of degree zero in z, z or there exists 7 and j,j € such 
that n-j — 1 and m~- = 1. In the first case one has that 

{ATeX K C L } - {Af E (z,z),u k u 1 } =0 
In the second case, then 

{Af E (z,z),u k u l z n z m } = u k u l ^El 1 ^ j {z 3 z ]1 u k u l z n z m } 

7 ieJ 7 

= u k u l z n z m El 1 J2 ~ mj) = 

□ 

Thus the _E-norm is invariant under the dynamics of a noncoupling Hamiltonian 
of degree at most 2 in z, z. In the more general case one has 

Corollary 6. Let Z be a non coupling polynomial. Assume that it has a smooth 
vector field, then there exists C such that 

\{M E ,Z}(u,u,z,z)\<CW E (z,z)f 2 (63) 

3.3. Normal form construction. Let us rewrite £ (ip, ^) as follows 

£ = Ho + eV, 

where 

Ho(V^):=«(ICi| 2 + IC 2 | 2 ) + E^IO| 2 , V:=±±^, (64) 

and 

VM,$) ■■= -(K2I 2 - ICil 2 ) + (65) 
e 

We are going to prove that there exists a canonical transformation T t which gives 
the Hamiltonian the form 

£ o % = Ho + eZ c + TZ (66) 
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where Z e is a non coupling polynomial and 1Z has a smooth vector which is expo- 
nentially small with e -1 . 

The construction will be recursive. To this end we assume one has been able to 
construct a canonical transformation T r putting the Hamiltonian in the form 

£o%. = E {r) =H + eZ< r > (e) + e r+1 U (r} (e) (67) 

with Z^ being a noncoupling polynomial and where TZ^ r \e) has a vector field 
which is bounded, uniformly with respect to e (Z^ — and 7?/°) = V e ). We 
look for an auxiliary Hamiltonian Q r +\ such that considering the corresponding 
Hamilton equations 



and the corresponding flow 
?' + 1 in place of r. 
Explicitly one has 



f/ r+1 one has that o (j) e r+l is in the form 1)67(1 with 



£ (r) o^l = Ho + eZ^ (68) 

+ e r+1 {n^ + {n ,g r+ i}) (69) 

+ Ho o <t>Cl -H - e r+1 {H 0i g r+ i} (70) 
+ e(zMo#£-ZM) (71) 

+ e r+1 (n^ o cjf^l - (72) 

Now, it is quite easy to understand that (|70H72(I are higher order terms (they will 
be estimated later), while 1(69(1 is the term of order e r+1 . Thus, if one is able to 
choose g r +i so that 

n {r) + {H ,g r +i} 

is non coupling, then the coupling terms are pushed to order r + 2, and 



(73) 



7 r+l 



is non coupling. 

The main step of the proof is the construction and the estimate of such a function 
Q fulfilling the homological equation associated to ((73(1 . 

3.4. Framework and notations. Before proceeding to the construction and the 
estimation of such a Q it is useful to extend the setting in which we will work. 
Indeed, since we are working with analytic functions it is useful to work in the 
complexification of our space. More precisely, we consider now a phase space in 
which the variables v :— u and w := z are independent of u, z. Actually this is 
equivalent to complexify the space in which vary the variables p, q of section 13.11 
With such an extension the Poisson bracket will take the form 



{H,IC} = 



2 

fc=i 



dU dJC dJC OH 



oo 
j=3 



dH dK. OK dU 



dzj dwj 



dzj dwj 



duk dvk duk dvk 
Such a phase space will be denoted by X£ . In Xq we will use the norm 

\{u,v,z,w)] := max{||(u,u)|| c2 , IK^w)^} 



(74) 
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with ||(2, w)\\ E defined by 

\\(z,w)\\%:=\\zf E + \\wf E , where \\ ■ \\ 2 E = M E (■ r) ■ (75) 

Sometimes it is also useful to use more compact notations for the phase space 
variables, thus we will also write £ :— (u,z), r] :— (v,w), £ :— {CiV)- 

Moreover, given R > we will denote by Br the ball of radius R and center 
in the phase space: 

B R =«e^ : m < R} ■ 

Given an analytic function hi with Hamiltonian vector Xu analytic as a map from 
Br to X£ we will denote 

\X h \r := sup |"*w(01 and Mr ■= SU P 1^(01 > (76) 
with norm defined by (|74ll . 

Remark 17. Bi/ Corollary^ our initial perturbation V t , defined by (|65[) . *s sucft 

\X Vt \ R <^R+^R^ +1 (77) 

3.5. Solution of the Homological equation. In this section we will construct 
and estimate the solution of the Homological equation associated to l|73|) : that is 
we look for Gr+i such that (|73[) is non coupling. 

Precisely, consider the Gauge invariant Hamiltonian function TL and decompose 
it (Remark j!5|l) as TL = T + Z with Z non coupling and T coupling, and consider 
the equation 

{7~to,G} = J~, T= coupling part of TL (78) 
then we are going to prove the following. 

Theorem 5. If TL has a Hamiltonian vector Xn analytic as a map from Br to 
Xg, then there exists a solution Q of l|78|l which is coupling and has a vector field 
Xg with the same analyticity properties; moreover there exists a positive a = a(h) 
such that the following ineguality holds 

\Xg\ R <^\X n \ R (79) 

moreover a(K) > C^h?! 2 

Before starting the proof we need some preparation. Since T is coupling then it 
is of degree at most 2 with respect to z and w. Therefore, we can decompose T as 
follows 

T :=T W + T 2Q + +T m (80) 

with 

T 10 = (F 10 {u,v);z) 1 T 01 = (F 01 {u,v);w) (81) 
T 2a = (F 20 (u,v)z;z), T 02 = (F a2 (u,v)w;w) (82) 

T 11 = (F 11 {u,v)z-w) (83) 
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and (.; .) be the scalar product of the real Hubert space I 2 defined by 

z^;z^):=^z^zf (84) 

3>3 

and F 01 ,F 10 which are X£ valued functions of (u,v) and F 20 , F 02 , F 11 functions 
taking values in suitable spaces of linear operators. We will denote by F? 1 ^, v) the 
components of F 01 ; we will denote by F® 2 (u,v) the matrix elements of F 02 (u,v), 
and similarly for the other quantities. 

Lemma 4. The Homological equation ()78(l has a formal solution Q with the same 
structure as T (cf. H8Q(I - (I83(I ), but with the quantities G's defined by 

F}°(u,v) m F^f^v) 

G ^ v ^-W^r g >' v): =^aH) (85) 

G » M := l (x J + x l -2ny G ^ u > v):= - i (x j + A l -2ti) (86) 

F n (u v) 
G)}{u,v):= J ' V \\ (87) 
i{Xj-Xi) 

Remark 18. In (jH3) one always has j ^ I since terms of the form F^(u,v) are 
always noncoupling; furthermore, one has that j and I cannot belong to the same 
set J 7 . The solution Q is coupling. 

Proof. We consider explicitly only the term Q m , since all the other terms can be 
studied exactly in the same way. First define 

H^'(u, v) := fl(uivi + u%v 2) , H^°°\z, w) :— XjZjWj, 

3>3 

so that H = H {1) + H {ca) . Now one has 

{n ,g 01 } = {n^,g 01 } + {n^,g 01 } (88) 



J2[{n^,Gf}-ix jG ' 



Wj 



where we simply computed explicitly {Ti^°°\ G 01 }- Thus to ask 

{Ho,G 01 } = F 01 (89) 

is equivalent to ask 

{H^,Gf]-iXjGf = Ff (90) 

We compute now the above Poisson Bracket. To this end decompose G" 1 in Taylor 
series, one has 

Gf{u,v)=Y. G tki^ 1 . (91) 

kl 

from which 

{H^,Gf} =J2G^m\ - \k\)u k v l (92) 

kl 
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but, assuming that Q 01 is coupling (which will be verified in a while), due to the 
limitation (|58|l . one has |fc| — \l\ = \m\ — \n\ = —1 (due to the fact that Q 01 is linear 
in id), and therefore 



01 



(93) 



Inserting this expression in l|90[) one gets 



■01 



iVLGf - iX 3 Gf = 

which shows that the function Q 01 of (|85|1 actually fulfills the equation (|89|) and is 
coupling. All the other terms can be studied in the same way; a detailed proof is 
omitted. □ 
End of the Proof of Theorem^ We will explicitly prove only the estimate of the 
norms of the vector fields of Xgm and Xgn , the other being similar and simpler. 
We start with Xgoi . Consider first the vector field of J- 01 , whose components have 
the form 



/dF 01 



\ dvi 



(94) 

X^ 01 = iF UL (95) 
Since T is the coupling part of TL then 

X^ox =d w X%(u,v,0,0)w (96) 
and therefore, adding also the v components, Cauchy inequality implies that 



dF 



01 



d(u, v) 



(u,v) 



1 



< -= sup[X^' >(u, V,Z, li- 
ft B R 



( """ ! , „.)- = L 

R 



X, 



(u,v) 

u 



(97) 



where the norm at l.h.s. is the norm as a linear operator from to C 4 . Using 
this inequality and making use of (|85|l it is very easy to estimate X^ '^: to this 
end define Wj := Wj/i(Xj — il) and let C > 1 be such that 



C^h := inf \X 3 

i>3 



01 



(98) 



Then one has 

\\w\\ E <C\\w\\ E /H 
and, using the definition of Q 01 , cf. H85|l . one has 



/dF 01 



and similarly for the v components, and thus 



Y (u,v) 



c 

< — 

r h 



x. 



{u,v) 
H 



(99) 



(100) 



for some C > 1. The estimate of the z components of the vector field is simpler 
(due to the simpler form of the z component) and collecting the two one gets 



C 

\Xgoi\ R < — \X H \ R 



(101) 
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We come now to the estimate of the vector field of Q 11 . Preliminary to this 
estimate we remark that the components of the vector field of J- 11 are given by 



= 4 w X%' v >(u,v,0,0)[z,w] (102) 

X™tt = d z X%(u, v, 0, 0)z = F 1Y z (103) 

thus in particular, by Cauchy inequality, 

supHF 11 ^,*;)!!^^ <hx%\ R (104) 



where the norm at l.h.s. is the norm as a linear operator from £ E to i E . 

In order to estimate XZu we have just to estimate the norm of the operator 
(from t E to itself) whose matrix is defined in l|87|) . To this end remark that the 
boundedness of G 11 as an operator from l 2 E to itself (or to its dual) is equivalent to 
the boundedness of the operator with matrix elements tiG^Sj as an operator from 
I 2 to itself, where ti, Sj are suitable positive numbers (in fact, £ 2 and £ E are spaces 
with not equivalent norms). Thus Lemma|S]of the appendix ensures that 

S up\\G n (u,v)\\ ee <-\X^\ R , a>C-^ 3 / 2 , (105) 

for some C > 1, since |Aj — A/ 1 > C _1 ft if j and I belong to different sets J 7 . 
Working in a similar way for the other components and the other parts of the 
vector field of the function Q one gets the result. □ 
From the proof (especially l|9t)|) . I|102|) and (|103fl ') also the following useful Lemma 
follows 

Lemma 5. The following estimates hold 

\Xr\ R <5\X H \ R , \X z \ R <6\X n \ R , \Z\ R <6\H\ R (106) 

3.6. Quantitative estimates. First we fix a positive R in such a way that X-p c 
is analytic on Br, and choose constants P, P* (that depend on h and on all the 
other parameters) such that 

\X Vr \ R <P, \V e \ R <P*. (107) 

All along this section we fix a small value of h, and will make explicit estimates 
so that at the end it will be possible to insert the dependence on h of the final 
estimate. 

Lemma 6. (Iterative Lemma) Consider a Gauge invariant Hamiltonian of the form 

£ {r] = Ho + eZ^ + e r+1 K {r \ (108) 

with non coupling and where = V t and Z^ = 0. Fix S < R/(r + 1), 
assume that the Hamiltonian vector fields of Z^ and oflZ^ are analytic on B^ r s, 



NORMAL FORMS AND NLS 



23 



id that 



X Z(-)\ R -r8 ^ 







if r = 
if r > 1 
if r = 



r(r) 



< 



R-r5 



< 



P 



J}* 

< — , 

R-rS eg 



with 



eo 



75P 



, where a = a(fr) > C 1 ft 3/2 , 



(109) 

(110) 
(111) 

(112) 



then, if — < 1/2 there exists a Hamiltonian function Q r +\ analytic on Bji-Sr 

generating the canonical transformation such that £^ o cff r+l has the form 

(I108|) and satisfies the estimates (|109D - llllfl with r + 1 in place of r, moreover the 
new Hamiltonian is Gauge invariant and the canonical transformation fulfills 



sup 
\£]<R-(r+i)S 



C-^" +1 (C) 



< e 



L " 

e J a 



(113) 



Proof. Decompose TZy) into its coupling part T r and its noncoupling part Z r . 
Define Z^ r+1 ^ := Z^ +e r Z r and use LemmaElto estimate \X zir+ i) \ R _ s , r+1 y Use 
Theorem |21 to construct Gr+i as the solution of 



then 



{Ti-O: Qr+l} — —Tr 



111, p 

\ X ^+i\R-rS ^ ~ \ X KM\R-r8 < ~ F 



(114) 



The Hamiltonian £^ o tjf was computed in subsection 13.21 and is given by equa- 
tions H68(l - I|72|l . which has the form 1)108(1 provided one defines 



n (r+l) = n (j-+l) + + n (r+l) 



where 



K^ +1) =e-^m, K ] = e- {r+2) m and 7^ +1 > = e^ 2 ^. 
Then, from Lemma fTTl with /i = e^ r+1 \ it follows that 



X n[^ 



(R-rS)-S Se r+ 



- J^+l^ \ X Z<~ r) \r-t5 \ X Q( r+ i 



) \R-rS 



< 



< 



r-1 



6^E 



p 

ac! 



60P 2 
aSe^ 
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Similarly 



X TZi r 



(R-rS)-S 



< 



< 



< 



M \ x km\r 



5c 
5£ 

5P 2 



2 r aSe r 



t5 

ae r 



(r+l) \n-rS 



< 



5P 2 
a5eg 



Furthermore, from Lemma ll2l with Tt = IZ^ and fj, 

25 

< - 

(R-rS)-S 



Ar+1) 



, it follows that 



X 



< 



< 



Se r+2 

25e r 

25P 2 
2 r aSe r n 



f\X 



KM \R-rS 
P 

ae r 



-rS \ X Q 



< 



(r+l) \R- r S 

25P 2 



Hence 



\ X n(r+v\fR-r5)-8 



< 



< 



P 

a 6 
P 
aS 



25P 60P 5P 



P 



2 r e r n 



^(60+^ 



(115) 



If r = then the second term is not present and the square bracket in (|115|l is 
equal to 30P; if r > 1, one has that the square bracket in (|115|l is less than 75-P/ep 
and thus the thesis on the vector fields follows. The estimates of the moduli are 
obtained in a similar way from Lemma [5j □ 
It is clear that the Hamiltonian Ij24(l of the NLS fulfills the assumptions of the 
Lemma with r = and thus the Lemma allows us to put our Hamiltonian in normal 
form up to any order r. To obtain the exponentially small estimate of the remainder 
take 6 = R/2r in order to fix the domain of definition of the final Hamiltonian, and 
then choose an optimal value of r, this can be done by minimizing the estimate of 
the remainder or simply choosing 

Ra 



150ePe 

(with [.] denoting the integer part). One thus obtains the following Theorem which 
is the main technical result of the paper, where Z — eZ^ r *^ and 1Z = e r * +1 TZ^ r *K 

Theorem 6. Consider the Hamiltonian £ (cfr. 124|) ). define 

aR 



150eP' 



(116) 



assume \e\ < e*/2 then there exists an analytic canonical transformation T such 
that 



£oT = H Q +Z+K, 



(117) 
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with Z non coupling; both Z and 1Z have a vector field which is analytic in -Br/2 
and fulfill the estimate 



\X z \ R/2 <12eP, |X TC | R/2 <ePexp[-(^) 



\Z\ R/2 < UeP* 



\n\ 



ft/2 



< eP* exp 



(118) 
(119) 



Moreover the transformed Hamiltonian is Gauge invariant and the canonical trans- 
formation fulfills 



sup R - T(01 < 2e- 

m<ft/2 a 



(120) 



It is also important to reformulate the Theorem computing P in terms of the 
original quantities u>, e, h i.e. using 1)77(1 . thus one easily gets the following 

Corollary 7. Fix a positive R, consider the Hamiltonian £ (cfr. 1(24(1 ). define the 
small parameter 



(121) 



i/ien i/iere exists > independent of H,u>, e, such that, if \x < /i*h 3 / 2 /2 i/ien i/iere 
exists an analytic canonical transformation T which transform £ into 11171) . where 

^3/2" 

|X z | H/2 < C/x , |X K | fl/2 < C^exp ' 



I ft/2 



< C/i, |^|r/2 < C/xexp 



(122) 
(123) 



Moreover the transformed Hamiltonian is Gauge invariant and the canonical trans- 
formation fulfills 



sup R - T(£)l < C- 

rei<ft/ 2 



(124) 



Theorem |3 is a direct corollary of Corollary [7] 

The manifold z = is approximately invariant for the dynamics of the system 
(|117|) and on such a manifold the dynamics is that of a Hamiltonian system with a 
Hamiltonian function which is an exponentially small perturbation of 

K.(u,u) := Z(u,u,Q,0). (125) 

Such a system has the additional integral of motion 

!(«,«) = £|«i| 9 . (126) 

3 

In the true nonlinear system one has 

Theorem 7. Assume fi < (i*h 3 / 2 and consider the Cauchy problem for the system 
((117(1 which is equivalent to NLS. Define 





r m^ 3/2 " 




max | exp 




.Nil'} 









(io 6em<7 i/ie initial datum for z) then, one has 

\\z{t)\\ s <C5 



(127) 
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\X{t)-X(0)\ <cs 2 

for all times t fulfilling 



\K(t) - K(0)\ <C[i6 2 



\t\ < 



(128) 



(129) 



Proof. First remark that due to the equivalence of the E norm and the s norm one 
has 



v^(Co) 

Thus, by ^ and one has 



dt 



<C f i(Af E ) 3/2 + C/iexp 



ii 



Me < Cn(Af E ) 



3/2 



which can be solved giving the estimate 

N E {t) < 2S 2 



(130) 



(131) 



(132) 



for the times (|129J) . The estimate of ||.z(i)|| s follows from equivalence of the norms. 
To obtain the estimate of the diffusion of X simply remark that M = X + Me is 
an exact integral of motion and use the estimate (|132|l . Finally, to estimate the 
diffusion of /C remark that the total Hamiltonian is an integral of motion, but one 
has 



\£oT-H + K,\ < Cfil \\z{t)\\' a + exp 



(133) 



which using the previous estimates implies the thesis. □ 
Proof of Theorem^ Define first M. :— T(<I>o). I n order to distinguish between 
the original coordinates of the NLS equation and the new coordinates introduced 
by the transformation T we will denote the new variables by adding a prime, i.e. 
we will write £ = T((^'). Remark that since T is Lipschitz together with its inverse, 
then for any couple of points £, £ one has 



c-c 



< 



c-c' 



< c 



c-c 



and therefore 



d(ip°,M) <6=>d(( ',$o) <C6 ^ \\z' Q \\ s <CS 



(134) 



(135) 



where by a slight abuse of notation we wrote ip° = X '(C')- Thus eq. i|127[l implies 
(|21|l for the considered times. To get (12211 . just remark that (again with a slight 
abuse of notation) 



and thus (|T2T|l implies (J23). □ 
The proof of Corollary is also a direct consequence of eqs. H128|l when one 
considers the deformation induced by the change of variables. A detailed proof is 
omitted. 
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Appendix A. Proof of Theorem Q 

For possible future reference, in this section we will work in R d . 
We start by recalling some notations and definitions from Robert 18; that will 
be used in the following 

Definition 4. A function m : M. 2d — > [0, +oo) is called a tempered weight if there 
exist Co, No > such that 

m(X) < Com{X l ){l + \X ~ X^ , VX,X 1 £R 2d (136) 

Definition 5. A function a £ C oc (R 2d ) is called a semiclassical symbol with weight 
(m,p), with p>0, if for any multi index a there exists a constant C a such that 

\d a a(X)\<C a {i ™^ )plal , VIet M (137) 

In this case we write a £ E™. 

In the following we will also need an extension to H dependent families of symbols. 

Definition 6. A smooth map h t— > a(K) £ E™ is called an h-admissible symbol of 
class (m,p), if for any integer N large enough one has 

a(H) — ao + Hax + h 2 a2 + ... + ft a^ + r/v+i (138) 

with aj S E™' -2 - 7 and rpf+x(fi) bounded in E™' 2 ( N+1 \ Here E™'^ is the class of 
symbols with weight 

(m(l + \X\) K »,p). 

Given an ^.-admissible symbol a of class (m, p), one defines the corresponding 
Weyl operator acting on L 2 by 

Op w (a)^(x) := f e lL ^a fel^, h) ^(y)d d yd d £ (139) 

By the theory of JS] one has that, for any H > 0, such an operator is well defined on 
the Schwartz space. Under suitable conditions it extends to a selfadjoint operator 
on L 2 . 

Definition 7. A strongly admissible operator of weight (m, p) is a C°° application 

A : (0,fr») -> £{S{R d ),L 2 (R d )) 

such that there exists an admissible symbol a £ E™ such that A(h) = Op w (a(h)) . 

One of the most important properties of strongly admissible operators is the 
given by the following Theorem. 

Theorem 8. (Theorem 11-32 of HB]J Let A = Op w (a) and B = Op w (b) be two 
strongly admissible operators of weights (to, p) and (n, p) respectively; then AB is 
a strongly admissible operator with symbol c(fi) and weight (mn,p). Moreover, if 
a ~ ^2j>o fc Ja j an d h ~ Sj>o h J bj then c ~ zCj>o ^ 1 c o w ^ 

c i = E rrgfla) U WfitfaOWfiW (140) 

\ a \ + \f}\+k+l= 3 \ / V / 
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(141) 



where D x — —id x . Furthermore, for any N > 0, one has 

N 
3=0 

and, for any a, (3 there exists a positive finite M such that the following estimate 
holds (for simplicity we restrict to the case p = 0) 

d"dlS N+1 (x, 1 < Cm(x, On(x, (142) 

N 

^ (PmM(aj)pnM(rN+l(b)) + PnM {bj )p m M (tn+1 (a) ) ) (143) 

j=o 



+ ^2 PmM(aj)PqM(h) J rP m M(r N+1 (a))p nM (r N+ i(b)) 

N<j+k<2N 



(144) 



where we denoted by rjy+i(a) the remainder of the asymptotic expansion of a trun- 
cated at order N and we denoted 



PmM{a) 



sup 



a\ + \f3\<M, (x,£) m \ x ,t,l 



(145) 



Fix a positive integer s > 1 denote A — Hq (where Hq is the operator J2Jl), from 
Theorem |H1 it follows that A = Op w (a) with a suitable a having principal symbol 
a = (£ 2 + V) s . Denote also b :=£ 2s + V s and B = Op w {b). Since V(x) > 1 then 
one has 



^ < ^b(x,0 < a (x,0 < Cb(x,Q, V(z,0 e R 2d . 



(146) 



Lemma 7. One has 

\\Ai/>\\ L , < C \\B^j\\ l2 , € D(B), (147) 
\\BiP\\ L2 <C\\AiP\\ L2 , V^efl(4 (148) 

Proof. Consider &~\ then b~ x E £jf\ Denote 6 := Op™^ 1 ). By Theorem 
one has 

AB = Op w (y) + hO P w (S 1 ) (149) 
BS = I + /iOp w (<52) (150) 

with #1,2 estimated by (|142fl - <|144[) . Since ^ is bounded together with its deriva- 
tives, it follows that Ai := Op w {8\) is bounded. Similarly A2 := Op w (S 2 ) is 
bounded. 

Thus, using Neumann formula one gets that the operator I + hA 2 is invertible 
provided h is small enough. So, from 1)150(1 one has 

{I + hA 2 )- 1 BB = l B- 1 = {I + HA 2 )- 1 B. (151) 

Finally one has 

||AV>|| i2 = \\ABB-^\\ L2 < \\AB\\c^,V) P'^L* ( 152 ) 



< C\\(I + hA 2 )- 1 B^\\ < C\\BiP\\ L2 



□ 
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Appendix B. Technical Lemmas 

We start by a Lemma which is needed for the estimate of the solution of the 
homological equation. In its statement we will denote by £ 2 the real Hilbert space 
of the sequences (zj)j>3 endowed by the scalar product 

(z;z'} :=^2 z 3 z j=J2 £ z i z J 

j>3 7 j>3, jGJ 7 

Lemma 8. Let F : i 2 — > I 2 be a bounded linear operator, assume that the corre- 
sponding matrix elements Fji are different from zero only if j G J 7 and I € Jy with 
7 ^ 7'. Define a new linear operator G with matrix 

G * = ijj&sj <153 > 

Then there exists a positive C such that the following estimate holds 

C 

C (e 2 ,e 2 ) - -^3/2 W F Wc(i> 2 ,i 2 ) ' ( 154 ) 

Proof. First we recall that |Aj — A;| > C~ 1 h if j E J 7 and I S J 7 ' for 7 ^ 7', and 
that #J 7 < C/ft. Fix Z € J 7 /. First remark that 

Ei^i^ii^iiW) (155) 

then, by Schwartz inequality 

/ \ 1/2 / x 1/2 

X>,i< few £ isrrjjp 

Fix Z € J 7 ' and estimate 

J J 77^7 J^J-y J 

V 7=7'±l 7=1 7>7'+2/ \je./ T 1 Jl / 

but, due to the choice of the numbers i? 7 one has 

r C" 1 ^ if 7 = y±i 

\\j-\i\>\ Ey - E n if 7<7'-2 (158) 
[£ 7 _i-£V if 7>7' + 2 

thus (|156[1 is estimated by 



T * J ) p + % (W + ) (15<J) 

Since the sums in i|159fl are convergent due to our choice of the sequence E~ on has 

x - 1 C_ 
4- |A,-A 3 f " ft 3 
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which gives 



£i^h< Aiin^) ( 16 °) 



j 

^From this one has 



Gz\\%=Y. 



l.Zl 



I 



= (ei^i) e(e \g S i\ N 2 ) 



CIIFII 2 



□ 

We report now some Lemmas from [3] which are needed for the proof of Lemma 
Here p will be small parameter (not the small parameter used in the main part 
of the text) that in applications will be replaced by e r with some r. 

Lemma 9. Let Q : B p — > C, p > be a function whose Hamiltonian vector field 
is analytic as map from B p — > X£; fix a positive 5 < p. Assume p\Xg\ p < 5 
and consider the flow </>' of the corresponding Hamiltonian vector field. Then, for 
\t\ < p., one has 

\tf-l\ p s<p\X g \ p . (161) 
Proof. It is just an application of the equality 

m - m = J* §(*)*» = f*x g (z(8))d8 . 

□ 

Lemma 10. Consider Q as above and let TL be an analytic function with vector 
field analytic in B p , fix < 5 < p assume p \Xg\ < 6/ 3, then, for \t\ < p, one has 



-Ho^\ p _ s - ^ ^ ^\^\ P J \^n\ p 
Proof. First remark that, since 0* is a canonical transformation one has 

Xno^iO = ^-*(^(0)^k(^(0) (162) 

from which 

To estimate the first term fix 8\ := 8/3] we have 

sup r# - V(0)-in < ™p r#- { (e)-ii 

rsi<p-35 1 r?i<p-2<5i 

<^ sup fa _t (0-£l <f \x g \ p . 

Going back to 5 and adding the trivial estimate of the second term one has the 
thesis. □ 
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Lemma 11. Let Q and H be as above, fix < S < p, and assume p \Xg\ p < 5/3, 
then, for \t\ < p one has 

5 

I^wo0*-w| p _ 5 < ^ \Xn\ p H \Xg\ p 

Proof. One has 

XnoAO - x n(0 = (d0-*(^(O) - I) Xni^iO) + (0) - *w(0] • 

The norm of the square bracket is easily estimated using Lagrange Theorem and the 
Cauchy inequality in order to bound dXn . The other term was already estimated 
in Lemma ITU1 so we have the thesis. □ 
Although Ho has unbounded vector field the vector field of Ho o ^ - Ho is 
bounded, more precisely we have 

Lemma 12. Let H, Ho, T and Q as in H2.Ht that is Q is the solution of the 
Homological equation (17811 ; denote by (p 1 the flow of the corresponding Hamiltonian 
vector field and 

£(0 := W (^(0) - «o(0 - M {H , Q} , 
then the vector field of I is analytic and, for any S < p, satisfies 

\X e \ p _ 5 <^ T \X g \ p \X H \ p . 

Proof. One has 

w (^(0) - Ho(0 = f ^W(0)d* = 

where we used the homological equation (|78|l to calculate {Ho, G}- Then one has 

m= I" [H^{S,))-Hi)]dt; 
Jo 

Using Lemma [5] and Lemma ITTI one gets the thesis. □ 
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